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The connections between the Random Phase Approximation (RPA) and Many-Body Per- 
turbation Theory (MBPT) and its all order generalisation, the Coupled- Cluster Theory 
(CCT) have been explored. Explicit expressions have been derived for the electric dipole 
amplitudes for allowed and forbidden transitions induced by the parity non-conserving 
neutral weak interaction. The Goldstone diagrams associated with the RPA terms in 
both cases are shown to arise in MBPT and CCT and the numerical verification of this 
relationship is made for the allowed electric dipole transitions. 
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1. INTRODUCTION 

The Random Phase Approximation (RPA) has been successfully used in calculating 
core polarisation/relaxation effects in a variety of properties ^. It is an approximate 
many-body theory that has been formulated in a number of different but equivalent 
ways ^. However, its connection with MBPT and CCT has not been explored in 
detail to the best of our knowledge. We focus on this particular point in this paper. 
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The outline of the paper is as follows. We first derive an effective operator in 
MBPT for electric dipole (El) transition amplitudes for allowed and forbidden 
transitions induced by the Parity Non Conserving (PNC) neutral current weak 
interaction From the general expression in MBPT, we consider diagrammatically 
all the RPA diagrams pertaining to zeroth and first order in the residual Coulomb 
interaction. This is followed by theoretical understanding of these effects starting 
from Hartree Fock (HF) equations to show how these effects can be represented as 
linear equations which in turn can be solved to self-consistency taking the residual 
Coulomb interaction to all orders. At the end, starting from the basic formalism 
for CC method, we compare the RPA diagrams both theoretically and numerically 
with the corresponding diagrams in CCT. 

2. General form of an effective operator in MBPT 

In perturbation theory, the functional space for the wave function is separated into 
two parts; a model space (P) and an orthogonal space (Q). The basic idea of such 
a division is to find an effective operator which act only within the limited model 
space but which generate the same result as do the original operators acting on 
the entire functional space. Here in the sections which follow, we derive an effective 
operator for allowed and PNC induced El transitions using perturbation theory. 



2.1. General form of an effective operator for allowed El 

transitions 

We start the derivation with the total Hamiltonian as 

H = Ho + Ves (1) 

where Hq and Ves are the imperturbed one-electron and perturbed two-electron 
operators. If |\l/o) is the atomic state function (ASF), then it satisfies the equation, 

Hol^o) = Eol^o)- (2) 

where Eo is the energy eigenvalue of the ASF. Here we are interested in allowed El 
dipole transitions between the eigen states of the atomic Hamiltonian H as given 

by 

El = (^'/jIDI*^) (3) 

where a and (3 denotes two different ASFs of different parity and D the dipole 
operator. Considering fl as the wave operator which upon acting on an unperturbed 
part generates the exact state, the observable El reduces to 

El = {^f^\Q"Dn'\^^°^) (4) 

Here, and denotes the initial and final atomic state functions. Once O is 
known. El can be computed. Starting from Bloch equation^ and considering only 
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terms of different orders in Coulomb interactions, we get El with -De// operator 
defined as 



oo 

^e//=E[^ei'"']'^E[^-)] 



m=0 



(5) 



n=0 



where O™'" refers to the wave operator with different orders of the residual Coulomb 

interaction denoted by m and n corresponding to the initial and final states. We 
consider -De// to have only connected diagrams. By putting n = and m = 0, it 
can be verified that it reduces to unperturbed contribution. Hence, a very general 
effective operator can be rewritten as 



D 



(n) 
eff 



m=0 



/{n — m) 1 



Hence El to any order n is given by 

isi(") = (*(°)|D5)|*(°)). 

For El^^^ with one order in Coulomb interaction, -D^V, reduces to 



Dl]) = DR'VesP + {R"VesP)''D 



(6) 



(7) 



(8) 



where R' and R" are the resolvent operators given by 

l^7)(^7l 



E„ 



R" 



7^M 

1*7) 



E2 



En — El 



(9) 



Here a, /? are in the model (P) and 7 in the orthogonal (Q) space with E defining 
their corresponding eigenvalues and M being the dimension of the model space 
considered. Considering the ASFs {a, /3) to be single determinant with single open 
valence shell we consider diagrams only of the form given in Fig. 1. The zeroth 



Fig. 1. Form of the diagram for El with v and v' denoting valence hnes. 



order El diagrams are represented in Fig. 2. Here diagrams (2b) and (2c) denote 
the direct and the exchange parts respectively and diagrams (2d) and (2e) their 
Hermitian conjugate parts. Higher order El diagrams can be obtained by taking 
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(n) 

different orders of Dpjj terms and contracting by Wick's theorem to obtain diagrams 
only of the form given in Fig. 1. In the next section, we follow similar derivation by 
taking PNC as an additional perturbation along with Coulomb operator. 




f 



(d) (e) 



Fig. 2. Diagrams representing _El(^^ contributions with diagrams (2b) and (2c) denoting the direct 
and exchange contributions with diagrams (2d) and (2e) with their Hermitian conjugate part. Here 
square sign represents Dipole and the dotted Unas denotes Coulomb operator. Fig (2a) corresponds 
to the Dirac-Fock part. 



2.2. General form of effective operator for PNC induced allowed 
El transitions 

We start the derivation with the total Hamiltonian as 

H = Ho+Ves+ HpNC (10) 

where Ves and Hp^c are two-electron and one-electron operators and Hp^c is 
expressed as 

HpNC = ^QwJ2^5p{re), (11) 
^ e 

with 

Qw = 2 [ZCip + NCin] ■ (12) 

Here Z and N are the number of protons and neutrons respectively and Cip and 
Cin are the vector (nucleon) - axial vector (electron) coupling coefficients whereas 
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Gf is the Fermi coupling constant and p(re) is the normahscd nucleon number 
density. The matrix element of Hpmc scales as Z^. Treating Hp^c as a first order 
perturbation, an ASF can be written as a state of mixed parity. 

1^-) = + |5'(corre)) (13) 

where I denotes the unperturbed part which is even/odd under parity and 
\'^{corTe)) denotes the correction due to the perturbation which is opposite in 
parity with respect to the unperturbed part. Due to the mixing of parity in the 
ASFs, one can expect a non zero electric dipole transition amplitude between states 
of same parity denoted by EIPNC as given by 

EIPNC = {^0\D\^oc) (14) 

where a and (3 denotes two different ASFs. Let us consider that we have only 
one parity (even+/odd—) in the model space. Considering Q as the wave operator 
which upon acting on an unperturbed part generates the exact state, the observable 
EIPNC reduces to 

EIPNC = on' \^+<^^^) (15) 

where the single and double prime on denotes the perturbation on the initial 
and final states. Once f2 is known, EIPNC can be computed. Starting from Bloch^ 
equation, a very general effective operator for various orders of residual Coulomb 
interaction with one order in Hp^c can be derived similar to the previous case. 
Defining fies with various orders of residual Coulomb interaction and Q.pnc with 
one order in Hpnc and different orders of residual Coulomb interaction, we can 
derive the general effective operator as 

n 

^5/ = E ["ei'") + n'^^Pc]'^}:-'^^ + ^Pnc"'^ (16) 
m=0 

where m and n refers to various orders of perturbation on initial and final states. 
For n = 1 with one order in iJpjvc and zero orders of residual Coulomb interaction, 
/Jg^j reduces to 

£>f/^ = DR'HpNcP + {R"HpNcP)''D 

where R and R' are defined as in the above case. By considering the ASFs to be 
single determinant with single open valence line and the fact that D and Hp^c 
are single particle; operators, we consider diagrams only of the form given by Fig.l. 
All the possible zeroth and first order diagrams with the corresponding expressions 
are given in reference'*. Out of that we are interested only in the RPA kind of 
diagrams which can be represented as linear equation starting from Hartree Fock 
(HF) equation and solved to all orders of residual Coulomb interactions. 
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3. Core Polarisation effects (RPA) in MBPT 

Core polarisation effects arise from the residual Coulomb interaction which is 
treated as a perturbation in MBPT. At each order of perturbation there is a single 
excitation from the core. The remaining correlation effects involving multiple ex- 
citations will not be discussed in this work. In RPA^'^^^ theory, the core electrons 
get perturbed in the presence of an oscillating electric field. By taking these orbital 
modifications into account in DF potential leads to coupled equations for the elec- 
tric dipole perturbed functions. In the subsections below, we show the equivalent 
terms/diagrams for the above effects in MBPT. 



3.1. RPA effects in MBPT for allowed El transitions 

The El transition amplitude has been given earlier in Eq.(3). This transition am- 
plitude in the zeroth and first order can be expressed as 

^l(o) ^ ^^W|£,|^(o)) (17) 

and 

= + (18) 

respectively. Using the general diagrammatic rules, all the possible diagrams of the 
form given in Fig. 1 can be obtained. This is given in Fig. 2. Converting the diagrams 
to expressions, we get 

Elf^ = {f\D\i) (19) 

and 

Here, we define a and q to be core and virtual orbitals with their corresponding 
single particle orbital energies denoted by e. Tilde on T4s refers to the inclusion of 
exchange terms. From the above expression, it is clear that due to residual Coulomb 
interaction at first order the core orbital 'a' is excited to a virtual orbitals defined as 
'g'. Hence the above expression denotes the first order RPA contribution in MBPT 
for allowed El transitions. Defining oj = ef — Ci, the higher order RPA diagrams 
can be obtained by solving the recursive relation as given by 

^(n) ^ '(T- {fq\Ves\ia)E^ag~^^ {fa\Ves\iq)E^a~^^ 

aq ^ aq ^ 

where n denotes the order of residual Coulomb interaction. Determining the above 
equation to self consistency first for core to virtual amplitudes and then using it for 
valence to virtual amplitudes is equivalent to taking all the RPA diagrams of the 
kind represented in Fig. 3 to all orders. 
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Fig. 3. The kind of RPA diagrams which are solved to all orders by representing the effects as 
linear equation and solving to self-consistency. Exchange and Hermitian conjugate diagrams are 
also included. 



3.2. RPA effects in MBPT for PNC induced El transitions 

The influence by an external oscillating electric fleld on the single particle orbitals 

of an atom can be obtained by solving DF equation in that field. By taking the 
matrix elements of th PNC operator between these perturbed states gives rise to 
the electric dipole transition amplitude which we are interested in. In the lowest 
order, we get 

ElPTVC^i) = (/^lifpjvcN) + {f\HpNc\^) (22) 
where we have defined 

i(../n = E¥^^ (23) 

Here e/ denotes the intermediate single particle energy. With additional mathemati- 
cal manipulations (as given in reference'') one can derive the all order RPA equation 
which in turn can also be represented as a linear equation and solved to all orders. 
The lowest order RPA contribution is shown in Fig. 4. Similarly by taking one 
order in residual Coulomb interaction, the second order electric dipole transition 
amplitude takes the form 

EIPNC^^-^ = T if\HpNc\p){pm) (24) 
^ (ei+u-ep) ^ ^ 



p 



^ y {f\HpNc\p){pa\Ve.s\tq){q\D\a) 
{ci + u) - ep){ea + - eg) 

_ \- {f\HpNc\p){pq\Ves\ia){a\D\q) 
^ {a + u) - ep){eg + oj - Ca) 

with 'p,q' and 'a' refers to virtual and core orbitals. The above equation describes 
the first order RPA effect as it involves the excitation of the core orbital 'a' to a 
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Ap 



hi 

(a) 



A, 



(b) 



Fig. 4. Diagrams representing EIPNC^^^ expression for RPA effects. The square and circle denotes 
dipole and PNC respectively. 



virtual orbital 'q' through the residual Coulomb interaction. The difference lies in 
the presence of the Hp^c perturbation acting either on the initial or final states. 
Comparing diagrammatically, the first and second terms of Eq. (24) is equivalent to 




(a) (b) (c) 

(Pseudo diagram) 

(because the sense of a and q are different (diagram witli local energy 
from the accepted convention) denominator) 



Fig. 5. Diagrams representing EIPNC^'^^ expression for RPA elfects with square and circle sign 
representing dipole and PNC perturbations. The dotted line denotes Coulomb interaction. 



the RPA diagram in Fig. (4a) and (5a) respectively. The third term is represented in 
two different ways as RPA diagrams in Figs. (5b) and (5c). In the first case, in Fig. 
(5b), the sense of the core and particle is different from the accepted convention with 
respect to D vertex and hence it is called a pseudo diagram. Whereas in the second 
case, by taking the negative sign inside the expression for E\PNC'^'^\ Fig. (5c) can 
be interpreted by local energy denominator for both D and Coulomb vertex. 

Comparing the above terms with all the first order MBPT terms as discussed 
in ^, we can find that the pseudo/local energy denominator diagram can be ob- 
tained by adding two MBPT diagrams as shown in Figs. (6b) and (6c) respectively. 
Similarly, the Hermitian conjugate pseudo diagram can also be obtained by adding 
the corresponding Hermitian conjugate MBPT diagrams. Fig. (6a) represents the 
normal RPA diagram as in the Fig. (5a). 
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(a) (b) (c) 



Fig. 6. MBPT diagrams corresponding to RPA. MBPT(l) corresponds to the normal RPA diagram. 
MBPT(2) and MBPT(3) are added to get the second term of the E1PNC^^> expression for RPA 
or in other words the pseudo/local energy denominator diagram. 



Similar to the allowed El transition, in the case of PNC, the influence of HpNc 
as a perturbation can be treated in the framework of RPA. In the next section, we 
show the basic formulation for CC method and show how RPA terms arises in the 
above mentioned method. 



4. RPA effects in Coupled Cluster (CC) method 

We have so far been discussing about perturbation theory in which some of the 
terms of the perturbation theory were grouped together and evaluated to all orders. 
This is similar in spirit to the Coupled Cluster method where a particular class of 
MBPT diagrams is calculated to all orders in the residual Coulomb interaction. 
The many-electron wave function in CC method is given by 

|*o) = e(^){e«}|$o), (25) 

where |<I>o) is the reference state and T and S arc the cluster operators which 
considers excitations from core and valence to virtual orbitals. 

The Hamiltonian in the presence of the PNC weak interaction is given by 

H = Ha + HpNc (26) 

where Ha and HpMC denotes the atomic Hamiltonian and the PNC perturbation 
with Gp showing the strength of the perturbation. Hence 

Here we denote the T and S with a superscript '0,1' as unperturbed and as PNC 
perturbed cluster amplitudes. Once the cluster amplitudes are computed, EIPNC 
can be obtained by taking one order in Hppfc perturbation in initial/final states ei- 
ther through T 01 S cluster amplitudes ^ connected by the D operator. Whereas for 
the allowed El calculations, the unperturbed cluster amplitudes in the initial/final 
states are connected by the D operator. In the sections below, the RPA diagrams 
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f 




(c) (d) 

Fig. 7. Coupled cluster counter part of fig 2. Here the dotted line is replaced by a solid line 
representing an all-order Coulomb interaction vertex. 

are compared both diagrammatically and numerically with CC for allowed El tran- 
sitions. 

4.1. Diagrammatic comparison of RPA diagrams corresponding to 
CC diagrams for allowed and PNC induced El transitions 

We now consider CC diagrams corresponding to the RPA effects. In CC theory for 
allowed and PNC perturbed dipole transitions, we consider D as an operator. The 
Coulomb operator acts as a perturbation for the allowed transitions, whereas for the 
PNC induced El transitions, we have PNC perturbation in addition to it. Therefore 
we need to first find out from which side of the dipole operator the perturbations 
act. 

The diagrams for the allowed El transition amplitude in the framework of CC 

theory has been discussed in one of our earlier papers ^. The RPA diagrams given 
in Figs. (2b-2e) correspond to the DS2 and S^D diagrams of CCT (See Fig 7). In 
the case of PNC induced El transitions, the MBPT diagram (4a) and (4b) in Fig. 
4 corresponding to zeroth order RPA effects can be got from '"^ where the 
superscripts denotes the order in PNC and the residual Coulomb interaction in the 
order of preference. By looking at the regular RPA diagram in Fig. (5a) contributing 
to first order in Coulomb and PNC, we find that it is equivalent to 082^'^^ dia- 
gram and its Hermitian conjugate in CC theory. Whereas the pseudo/local energy 
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denominator diagram, which is shown as equivalent to addition of MBPT diagrams 
in Fig. (6b) and (6c) can be got from the term 5j^^''^^£'S'2^'°^ and its Hermitian 
conjugate respectively. 

By comparing the terms in CC method for allowed and PNC induced El tran- 
sitions, we find that the additional diagrams arises solely due to the presence of 
PNC operator acting either on the incoming or the final vertex. We find that PNC 
acting on the same vertex as Coulomb operator, leads to the regular diagram and 
PNC and the residual Coulomb interaction on the initial and final vertices leads to 
the additional pseudo/local energy diagram. 

5. Numerical comparison of RPA diagrams corresponding to CC 
diagrams for allowed El transitions 

In order to demonstrate that the RPA effects are contained in CCT we have com- 
pared the all order RPA results with the sum of DS2 and 5*2-0 contributions in 
our CC calculations for the allowed 3s^6'i/2 — > 3p^Pi/2 and 3s^5i/2 — > ^P^P3/2 
transitions in Na^. We first generate the single particle DF orbitals with Na"'"(2p^) 
as the starting potential. The GTOs were generated using the Finite Basis Set Ex- 
pansion (FBSE) method with a primitive basis set consisting of 35s-32p-25d-25f. 
With appropriate energy cutoffs in the discrete core and continuum virtual orbital 
spectrum, the calculation was done with lls-10p-9d-8f basis. 



3s - 




3s - 




RPA 


RPA-CC 


RPA 


RPA-CC 


0.044 


0.044 


-0.061 


-0.062 



Table 1. Comparative results of El reduced matrix elements (in a.u.) arising from RPA and 
RPA-CC {DS2 + S^D) for 3s — > 3pi/2 and 3s — > 3p3/2 transitions in Na+. 

As stated earlier DS2 and terms in CCT (Fig. 7) correspond to the RPA 
diagrams given in Figs. (2b-2e). In table 1 we give the numerical values of those two 
contributions for 3s — > 3pi/2 and 3s — > 3p3/2 transitions in Na+ . As expected 
the agreement between the two results is very good. 

6. Conclusion 

Using analytical and diagrammatic techniques we have identified the RPA effects 
in allowed and parity non-conserving electric dipole transition amplitudes. We have 
demonstrated that these effects arise in MBPT and CCT. Indeed the all order RPA 
contribution is subsumed in CCT. 
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